Spontaneous spin ordering of Dirac spin liquid in a magnetic field 
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The Dirac spin liquid was proposed to be the ground state of the spin-1/2 Kagome antiferromag- 
nets. In a magnetic field B, we show that the state with Fermi pocket is unstable to the Landau 
level (LL) state. The LL state breaks the spin rotation around the axis of the magnetic field. We 
find that the LL state has an in-plane 120° q = magnetization M which scales with the external 
field M ~ B a , where a is an intrinsic calculable universal number of the Dirac spin liquid. We 
discuss the related experimental implications which can be used to detect the possible Dirac spin 
liquid phase in Herbertsmithite ZnCu3(OH)6Cl2- 

PACS numbers: 75.10.Jm, 75.50.Ee 



Spin liquids (SL) , defined as the ground states of spin 
systems with half integer spins per unit cell which does 
not order magnetically and/or break translation symme- 
try, are believed to contain fundamentally new physics 
beyond the Landau's symmetry breaking characteriza- 
tion of phases. After years of search, a promising can- 
didate finally emerged in the spin-1/2 Kagome system 
Herbertsmithite ZnCu 3 (OH) 6 Cl 2 0,lEl3. Despite an an- 
tiferromagnetic exchange J sa 170 ~ 190KJ4L [5|, the sys- 
tem does not order down to 50mK[l|, 0, S|. Theoreti- 
cally the Dirac-SL is proposed0 as the ground state of 
the nearest neighbor Heisenberg antiferromagnetic model 
on Kagome lattice. Unfortunately the spin susceptibility 
Xs is consistent with ~ 4% magnetic impurities [J], and 
possible Dzyaloshinskii-Moriya interaction has also been 
proposed[5[. These may explain the hump of the specific 
heat C around 2K and also obscure the C oc T 2 and 
Xs oc T behaviors predicted by the Dirac-SL. In this pa- 
per we propose a unique signature of the Dirac-SL in an 
external magnetic field, which can be used to detect the 
Dirac-SL phase in experiment and numerical simulations. 

The Dirac-SL experiences no orbital effect in the exter- 
nal magnetic field because it is an insulator. The Zccman 
coupling g^BB-J2i Si polarizes spin along the field direc- 
tion, breaks time reversal and breaks SU(2) spin rotation 
down to U(l). Let us denote the direction of B as z- 
direction (note that the xy-plane does not have to be the 
plane of the two-dimensional Kagome system). We find 
that the Dirac-SL will spontaneously break the remain- 
ing S z -U(l) and form a staggered magnetization M (see 
FiglUb)) in the xy-plane which scales with B: M ~ B a . 
The positive exponent a is an intrinsic universal number 
of the Dirac-SL phase and in principle calculable. 

This unique signature of the Dirac-SL can be compared 
with a regular co-planer antiferromagnetic (AF) ordered 
phase. In a small magnetic field B, the magnetization 
of the regular AF phase would rotate into the xy-plane 
to maximize the susceptibility along the B direction. As 
B is tuned to zero, the in-plane magnetization remains 
finite. In contrast, for the Dirac-SL we predict that M 



vanishes as B a . This suggests that the Dirac-SL can be 
viewed as an AF phase whose long-range AF order has 
been destroyed by the quantum fluctuations. 

We begin with writing down the the low energy theory 
of the Dirac-SL[6], which includes four flavors of fermions 
coupled with compact U(l) gauge field in 2+1 dimension 
(QED 3 ): 

S = j dx 3 [\{e\ tiu d l _ l a u ) 2 + ^ ^ +cr (d^ - ia^) r^V+o- 
+ 5Z^-' T ( 9 » ~ ia ^) T ^-<y] + ■■ ■ . (!) 

a 

where the two-component fermionic Dirac spinon fields 
are denoted by V'ic where ± label the two inequivalent 
nodes and a the up/down spins. 

The spinon carries spin-1/2 and charge-0, and couples 
to the external magnetic field only by the Zeeman effect. 
In the presence of a magnetic field, the simplest guess is 
that the Dirac points will change into Fermi pockets due 
to Zeeman splitting. We will call this state Fermi pocket 
state (FP). But, can there be other states whose energy 
may be lower? 

The FP state has many gapless spinon excitations near 
the Fermi pockets and in general not energetically favor- 
able. One natural way to gap out the Fermi pocket state 
is to induce (internal) gauge fluxes and develop Landau 
levels. Here we take advantage of the appearance of zero 
energy Landau levels when Dirac particles are subject to 
a magnetic field, as is well known in the recent studies 
of graphene0, Q. If the gauge flux is adjusted in such a 
way that the zero-energy Landau level is fully filled for 
up-spin, and fully empty for down-spin, then the spinons 
are fully gapped (FigQ}. We will call this state Landau 
level state(LL). 

In the following we show that LL state has lower energy 
than FP state. We first compare the energies of LL state 
and FP state at the mean-field level (ignoring the gauge 
fluctuation) with fixed S z polarization. We set aside the 
Zeeman energy which is common between the two states. 
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FIG. 1: The The spinons is gap less in FP state with an 
electron-like Fermi pocket for spin up spinons and hole-like 
pocket for spin down spinons. The spinons are fully gapped 
in the LL state. Each Landau level is doubly degenerate due 
to two inequivalent Dirac nodes. 

The density of spin imbalance is: 

An=^ = ^^ = 4.J-f^V, (2) 
A A 4vr V v F J 

and the mean-field energy density is found to be 

AeS = ^ = ^(An)l, (3) 

where iVj and iVj, are the number of up and down spins, 
AE is the energy increase compared to Dirac point state, 
A is the system area, B is the magnetic field, and vp is 
the mean- field Fermi velocity. The factor 4 in Eq.@ is 
from the fact that there are two spin degeneracy and two 
nodal degeneracy. We choose our units such that H = 1. 

To make sure the LL state have the same spin imbal- 
ance as the FP state, each Landau level should contain 
AnA/2 states, i.e. the induced gauge magnetic field b 
satisfies An = 2 • j-, where the 27r is gauge flux quan- 
tum. Since at the mean-field level the internal gauge field 
costs no energy, the energy of the mean-field LL state is 
just the difference between the sum of the energies of all 
negative landau levels, and the filled Fermi sea: 

= C(f)ft 3/ V = C(f)An 3 /V 
V2n 2 ' 

The same result was obtained in Ref 0, U3]. From 

Eq.dSHl) we have = = 1.134, i.e., the LL 

state has lower mean field energy than the FP state. 

We may also ask whether in graphene a spontaneous 
magnetic field may be generated by this mechanism, if 
electrons or holes are introduced by gating. The answer 
is negative and the difference is that the physical elec- 
tromagnetic field has an energy cost of / 2 . Furthermore 
the field is three dimensional and the magnetic field must 
form closed lines. The magnetic flux would be divided 
into domains of size R with opposite signs. The size R 
can be estimated by minimizing the total energy, which 
is the sum of the magnetic field energy ~ B 2 R 3 and the 
Landau level energy gain ~ —VFy/heB/c £>i? 2 /$o where 
$o = — is the flux quantum. The optimal R* is found 




FIG. 2: (color online) (a) We convert the data in Table U 
into energy density Ae and spin imbalance density An (units 
defined in text) , and plot them together with the fitted scaling 
form Ae = aAn 3/2 + bAn 2 . (b): The XY spin order in the 
LL state. We also present the relative phase of Sf operator 
computed by projected landau level state in unit of n (see 
Eq.©). Note that we fix the left-bottom corner site as the 
site i in Eq.©, and scan the site j for the nine sites depicted. 

to satisfy ~ (a^) 2 , where ot = j- = ^ s the fine 
structure constant. This means that the flux through the 
entire domain is ~ 10~ 7< J>o , which is not self consistent 
with the condition that there is at least one flux quantum 
through the domain to support Landau levels. 

To include the effect of gauge fluctuations and to go 
beyond the mean field theory, we will use the Gutzwillcr 
projected wavefunction to calculate the energies of LL 
and FP states. To obtain the Gutzwillcr projected 
wavefunction we first write down the mean-field Hamil- 
tonian on lattice: H mean = J2ij XijfLfj*' where f la 
are the fermionic spinons. The mean- field ground state 
I'S meanixij)) "with mean-field parameters Xij is a spin 
singlet. The projected wavefunction \^fprj{Xij)) = 
mean(Xij)) removes the unphysical states and be- 
comes a spin state; here P = Ili(l~ n iT n il) ^ s the projec- 
tion operator ensuring one fermion per site. The physical 
observables can be measured on \^ pr j(Xij)) by a Monte 
Carlo approach [ll|. The Dirac-SL is characterized by Xij 
such that \xij\ = X is a bond independent constant and 
there are it fluxes through the hexagonal plaquettes and 
fluxes through the triangular plaquettes [6[. Note that 
the projected Dirac-SL state has no tunable parameter 
since \ on ly gives the wavefunction an overall factor. On 
a 16x16 unit cell lattice, the energetics of the FP and LL 
states are given in Table HI 

Is the Ae ~ An 3//2 law still valid after projection? The 
answer is positive because the S z and energy are con- 
served and have no anomalous dimension. The energies 
in Table [Jean be thus fitted (FigJ^a)) as 

Ae££ = 0.33(2)An 3 / 2 + 0.00(4)An 2 , (5) 
Ae P L r l = 0.223(6)An 3 / 2 + 0.03(l)An 2 , (6) 

where the units are chosen such that the unit cell spacing 
a = 1 and J = 1 and we include the first order correction 
to scaling, the An 2 term. From the coefficient 0.33(2) in 
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AN 





8 


16 


rpFP 
^Pr.i 


-0.42865(2) 


-0.42798(2) 


-0.42688(2) 




-0.42865(2) 


-0.42817(2) 


-0.42732(2) 


AN 


32 


40 


56 


„FP 
_C/p rj 


-0.42327(2) 


-0.42131(2) 


-0.41638(2) 




-0.42493(2) 


-0.42346(2) 


-0.41994(2) 



TABLE I: The per-site energies of the projected FP and LL 
states in unit of Heisenberg coupling J on a 16x16 unit cell 
lattice with periodic-antiperiodic boundary conditions. AN 
is the spin imbalance. The FP states has no total gauge flux 
through the sample, while the LL state has AK fl ux quanta 
through the sample which is uniformly distributed in the tri- 
angular and hexagonal plaquettes. 



Eq. ([5]) we are able to read off the effective fermi velocity 
v* F = 3 '° 2 ^ 2) = 0.28(2) ^ by fitting to the free fermion 
result Eq.((3]). v F is almost twice the mean- field value 
Vf = vis with x = °- 221J P^I- For Herbertsmithite this 
means v* F — 4.9-10 3 m/s assuming J = 200K. This is very 
close to the fermi velocity we found by a projected band 
structure study [l3} , where only one particle-hole excita- 
tion is considered. The closeness of the two results imply 
that the gauge interactions between many particle-hole 
excitations may only give small corrections for energetics. 

From Eq. (|5|6|) we see that the LL state has a lower 
energy than the FP state. Therefore the LL state may 
be the true ground state in the presence of a magnetic 
field. Due the presence of internal gauge flux, one can 
easily see that the LL state breaks parity (mirror reflec- 
tion). Surprisingly, it turns out that the LL state breaks 
the S z spin rotation as well. To understand this let us 
consider the low energy collective excitations in the LL 
state. Since the spinons ip are gapped due to the Landau 
levels, the gauge field a M is the only low energy excita- 
tion. We note that the total Hall conductance of the 
spin-up spinons is +1 and the total Hall conductance for 
the spin-down spinons is —1. Thus the total Hall conduc- 
tance is zero and there is no Chern-Simons term for the 
gauge field a^. As a result, the dynamics of the gauge 
field is controlled by the Maxwell term f^ u , which arise 
from integrating out high-energy fermions, and gives rise 
to a linear dispersed gapless photon mode. 

To understand the meaning of the gapless mode, we 
note that due to the non-zero Hall conductances ±1 for 
the spin-up and spin-down spinons, a flux quantum of the 
gauge field will induce one spin-up spinon and minus 
one spin-down spinon. Thus a flux quantum carries a spin 
quantum number S z = 1. The "magnetic field" of is 
nothing but the S z density. Therefore, the linear gapless 
mode of a M is actually the density fluctuations of S z . 
The appearance of a linear gapless mode of a conserved 
density as the only low lying excitation implies that the 
corresponding symmetry is spontaneously broken in the 
ground state |l4t HBj Thus the LL state contains an XY 



order: (S + ) ^ and the gapless gauge field is the 
Goldstone boson mode of the U(l) symmetry breaking. 
The "electric charge" of gauge field corresponds to 
the vortex in the XY ordered phase. This duality 
17J between the spin degrees of freedom and the gauge 
degrees of freedom allows us to translate the physics in 
the two languages back and forth. 

For example let us consider the monopole insertion 
operator V' which inserts a 2n gauge flux. V' flips 
one spin as mentioned before and can be expanded as 
V' ~ Yli e Wi Sf , where the summation is restricted to 
the area where the 2-k flux is inserted. Because 
has long-range correlation in the Coulomb phase, by the 
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above duality, the spin correlation 
also long ranged. 

To study what kind of the XY spin order that the 
LL state has, we have calculated spin-spin correlation 
function (S x ^S x j + S yi iS y j) of the projected LL state 
and found the XY ordered pattern as shown in Figl^b). 
This order pattern is referred as the q = magnetic order 
in literatures [ijj]. To confirm this is indeed the order 
pattern we directly compute the relative phases of Sf~ 
operators in the monopole insertion operator V'; 
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Wn 



1 flux quantalPS^Pln flux quanta) \ 
1 flux quanta|PS' l + P|n flux quanta) / ' 

(7) 



where \n flux quanta) is a LL state with n flux quanta 
and P is the projection operator. To understand this 
formula one should note that P\n + 1 flux quanta) = 
WP|n flux quanta), where inserts an extra 2-7T flux 
uniformly while flipping one spin, and thus ~ 
'^2 i e l9i Sf with summation over all sites. These relative 
phases can be computed analytically by choosing proper 
gauge (i.e., gauge that respects physical symmetries). To 
get the answer quickly we perform the Monte Carlo cal- 
culation on a finite lattice. For a 6 by 6 unit cell sample 
on torus with n = 4, we compute this relative phase for 
9 sites within the three adjacent unit cells and the result 
is given in Fig. [DJb). We find the relative phases are 
0, ±-j^ alternating between the three sublattices and this 
is exactly the 120° q = order pattern. [irjj] 

Because is in fact the magnetization operator, we 
conclude that the in-plane staggered magnetization M ~ 
B a scales with the external magnetic field, where the ex- 
ponent a is the scaling dimension of the monopole opera- 
tor at the Dirac-SL fixed point, because B is dimension-1. 
This exponent is computable by numerical simulations of 
QED3 or the field theory techniques such as 1 /N expan- 
sions, and subject to further study. 

So far we show that the FP state is unstable towards 
the Sz-broken LL state by numerical arguments. In the 
following we present an analytical argument that the S z 
symmetry is broken by studying the low energy effective 
theory. The FP state is characterized by a electron-like 
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Fermi pocket of spin up spinons fk\ and a hole- like Fermi 
pocket of spin down spinons fki . After doing a particle- 
hole transformation on the spin down Fermi pocket only 
fki ~ > h-*T we have and ft._fef carry same spin but 
opposite gauge charges. The "Coulomb" attraction be- 
tween the two particles will cause a pairing instability 
(/fct^-fet) = (fir fki) 0- This is a triplet excitonic 
insulator and this has been discussed in the context of 
grapheneflil liol]. Besides Coulomb attraction, there is 
also an Amperean attraction [2lJ between the currents of 
/fcl and h_ k ^ excitations. The "Coulomb" attraction and 
the Amperean attraction are cooperating for the same 
condensation. Note that fl^fki is a S z spin-1 object. 

Finally we discuss the consequences of this sponta- 
neous spin ordering in experiment. First the in-plane 
q = magnetization pattern and its scaling law M ~ B a 
is observable by neutron scattering. Second since the 
ground state breaks the parity and S z rotation symmetry, 
it is separated from the high temperature paramagnetic 
phase by at least one finite temperature phase transition. 
This transition can be first order or continuous. If the 
transition is a continuous one that restores the S z -U(l) 
symmetry, we expect it to be the Kosterlitz-Thouless uni- 
versality. Here we simply estimate the transition temper- 
ature T c ~ when [i B B < x Prj ■ This is because the 
external magnetic field B is the only energy scale if it is 
much smaller than the spinon band width x Pr ' 7 - Since 
these are intrinsic properties of the Dirac-SL in magnetic 
field, they can be used to experimentally detect the possi- 
ble Dirac-SL ground state in Herbertsmithite and other 
materials where a Dirac-SL may be realized. This XY 
spin ordering can also serve a way to detect Dirac-SL 
in numerical studies of the Kagome lattice Heisenberg 
model such as exact diagonalization. 

There are other proposals for the possible non- 
magnetic Valence Bond Solid (VBS) ground states in the 
spin- 1/2 Heisenberg model on Kagome lattice 22j, |23j. If 
one of those VBS states is the ground state of Herbert- 
smithite, it will surely not break S z symmetry in a small 
magnetic field because the VBS phases are fully gapped. 
The S z symmetry can be broken at an external magnetic 
field larger than the spin gap due to the triplon conden- 
sation. However the XY magnetic order generated in this 
fashion is unlikely to be the q = pattern because the 
VBS orders itself breaks translation and we expect their 
XY orders also breaks translation and have a large unit 
cell. Therefore the spontaneous spin ordering we present 
in this work can be used to differentiate the J7(l)-Dirac 
SL from VBS states experimentally. 

It is likely that in the current Herbertsmithite com- 
pound there is an energy scale below which the signif- 
icant amount of impurities and/or the Dzyaloshinskii- 
Moriya(DM) interaction start to play an important role. 
Recent experiments estimate the strength of the DM in- 
teraction to be ~ 15K[24J. In a strong magnetic field 
~ 30 tesla we expect that it is possible to suppress their 



effects and reveal the intrinsic property of the Dirac-SL. 
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We actually find that this relative phase is independent 
of n. Even when n — 0, namely the pure !7(1)-Dirac 
SL state, the relative phase is the same as the n — 4 
data that we presented up to statistical error. This indi- 
cates that the relative phase which we computed is an 
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intrinsic property of the £/(l)-Dirac SL, i.e., they are 
the spin-1 monopole quantum numbers. The Eq.([7]) is 
a systematic way to compute monopole quantum num- 
bers in any £/(l)-Dirac SL. A more detailed discussion on 
this method is presented as an appendix in Ref . [25(] . The 



monopole quantum numbers obtained in this method (see 
Pig. Efb)) are consistent with the results obtained by an 
independent study on the Dirac SL in zero magnetic field 
based on symmetry group analysis [l^]. 



